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I.  INTRODUCTION 

The  flight  stability  of  liquid-filled,  spin-stabilized  projectiles  heis  been  considered  for 
a  wide  variety  of  conditions.  Originally,  theories  and  experiments  were  centered  about 
the  case  of  large  Reynolds  number  {Re  =  a^0/i/).'^  ^  Many  complicated  non-steady 
effects  must  be  considered  for  practical  applications  if  the  Reynolds  number  is  large,  for 
example,  spin-up.^'  ■*  ’  ®  However,  if  the  Reynolds  number  is  very  low  {Re  <  50), 

the  effects  of  these  unsteady  processes  may  be  neglected.  It  is  then  possible  to  employ 
an  incompressible,  fully  viscous  finite-difference  solution  to  the  Navier-Stokes  equations 
for  a  fixed  precession  angle  and  steady  rates  of  spin  and  precession.  Pressures  hp  ’e  been 
measured  under  these  conditions  zmd  can  be  used  to  validate  numerical  simulations, 

This  report  reviews  available  methods  for  calculating  the  behavior  of  low  Re  payloads. 
Detailed  comparisons  between  two  finite-difference  simulations  are  made.  Conclusions  and 
recommendations  as  to  the  accuracy  and  computing  time  for  all  methods  are  presented. 


11.  BACKGROUND 

Two  such  stead}'  state,  incompressible  Navier-Stokes  (SS-INS)  codes  have  been  devel¬ 
oped.  Initially,  \  aughn  and  co-workers®  developed  an  explicit  code  for  the  steady-state  so¬ 
lution  of  the  dimensional  Navier-Stokes  equations  using  an  iterative  finite-difference  proce¬ 
dure  and  Chorin's  method  of  artificial  compressibility.  This  code  (SAND)  employs  uniform 
numerical  grids  in  the  radial,  axial,  and  azimuthal  directions  and  central  finite-differences. 
Recently,  Strikwerda  and  co-workers®  have  developed  an  implicit  code  for  the  steady  state 
solution  of  the  non-dimensional  Navier-Stokes  equations  using  an  iterative  finite-difference 
method  based  on  modified  line  successive-over-relzixation  (LSOR)  and  a  pressure  update 
from  the  gradient  of  the  velocity  field.  This  code  (UWISC)  employs  non-uniform  grids  to 
better  resolve  the  velocity  and  pressure  near  the  cylinder  wedls,  central  finite-differences  in 
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Journal  0/  Caiianrc,  Control,  and  Dynarrice,  Vot.  So.  3,  pp. 354-359,  May- Jane  1985. 
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BaUiitie  Research  Laboratory,  Aberdeen  Prating  Ground,  Maryland,  Deee'nber  1984.  (AD  150180) 
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Reieareh  Laboratory,  Aberdeen  Prating  Ground,  .Maryland.  February  1984.  (AD  139138)  AUo  AIAA  Paper  83-tl43,  Auguit 
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BaltnUe  Reieareh  laboratory,  Aberdeen  Proving  Ground,  Maryland.  BRL  .Wtmoraniium  Report  in  rrenaradon. 

Vaughn  H.  R.,  OberkampJ.  H’ .  and  Wolje.  tV.  R..  "Fluid  Motion  Iniide  a  Spinning  .Sutalin 
JcMrnal  oJ  Fl^td  Mechanics  Vpi.  ISO.  fp.  J9»S.  AI$o 


C'^itnder,"  Sarid\a  Hrpori  SAXD  Decemhtr  J96S. 


Spinning  Sutating  Cylinder 
AtimeWcaf  5o/uI»on  Jot  a  Spiti'ntng,  .Vufalin^  Fluid-Filled 


’Sinkuerda  J.  C.,  and  Sagei.  )  .  M..  ".•(  Sumerical  Study  0)  Flow  in  Spinniny  and  Conmg  Cylinder!,"  CRDC.SP-86007, 
Proceeding!  0]  the  tio3  Scientific  Conference  on  Chemical  De/enie  Reieareh.  Aberdeen  Proving  Ground,  .Maryland,  April 

1566. 


'•0 


the  radial  and  axial  directions  and  pseudo-spectral  differencing  to  represent  the  azimuthal 
dependence.  These  codes  have  been  used  in  the  present  effort  to  obtain  solutions  for  a 
wide  rzuige  of  dimensional  and  non-dimensional  parameters. 

In  addition  to  finite- difference  solutions  of  the  incompressible  Navier-Stokes  equations, 
analytical  solutions  to  the  linearized  equations  are  available.  Herbert^^  has  developed  a 
simple  model  of  the  viscous  flow  in  a  cylinder  of  infinite  length  that  is  spinning  and  nutating 
at  a  small  yaw  angle.  Due  to  the  infinite  cylinder  eissumption  of  Herbert,  the  theory  h^ls 
been  found  to  agree  with  UWISC  in  liquid  roll  moment  coefficient  for  aspect  ratios  in 
excess  of  3  (Figure  1).  For  the  liquid  side  moment,  coefficient,  Herbert’s  assumption  of  zero 
pressure  force  in  the  linearized  equations  has  yielded  inaccurate  results.  However,  since 
a  relationship  between  liquid  roll  and  side  moment  coefficients  has  been  established." 
Herbert’s  theory  cem  be  used  as  a  quick  and  effective  tool  for  initial  estimates.  Recently. 
Herbert  has  developed  a  spectral  collocation  method*^  for  the  Navier-Stokes  equations  and 
a  finite  length  cylinder.  Results  for  liquid  yaw  (side)  moment  are  in  good  agreement  with 
the  UWISC  code  (see  Figure  12  of  Reference  12). 

A  spatial  eigenvalue  method  has  recently  been  developed  by  Hall.  Sedney,  and  Gerber." 
The  Navier-Stokes  equations  are  written  in  an  inertial  reference  fraine  and  reduced  to  a 
set  of  linear  partial  differential  equations.  The  angle  of  coning  motion  is  assumed  small 
and  only  linear  departures  from  solid  body  rotation  are  considered.  To  obtain  boundary 
conditions,  the  motion  of  the  cylinder  walls  is  imposed  from  the  pioiectile  motion  which  is 
proportional  to  where  t  is  time  and  /  is  r,  the  non-dimensional  coning  frequency  for 

pure  coning  motion.  The  flow  variables  (perturbed  velocity  components  and  pressure)  are 
then  proportional  to  A  particular  solution  was  employed  which  satisfies  axial  and 

lateral  wall  boundary  conditions  but  not  endwall  conditions.  The  eigenN-alue  problem  is  de¬ 
fined  using  a  separation  of  variables  technique  from  which  an  infinite  sequence  of  complex 
eigenvalues  is  generated.  The  eigenvalues  are  determined  by  an  iterative  process  for  which 
sufficiently  accurate  initial  estimates  are  required  for  convergence.  The  flow  variables  are 
expressed  as  eigenfunction  expansions  whth  the  coefficients  determined  by  satisfying  the 
endwall  boundary  conditions;  a  least  squaires  and  collocation  method  have  been  used  for 
this  purpose. 


Comparisons  of  measured  liquid  moment  coefficients  with  spatial  eigenvalue  and 
UWISC  results  have  shown  the  consistency  of  both  methods.  How-ever.  since  spatial 
eigenvalue  methods  yield  results  in  significantly  less  computer  run  time,  they  are  per¬ 
haps  the  preferred  scheme.  Figure  2  shows  a  comparison  of  endwall  pressure  coefficient 
from  U\^  ISC.  S.4ND  and  the  spatial  eigen\’alue  method  with  experimental  values.'  Con¬ 
siderable  discrepancy  is  noted  for  the  SAND  code  while  the  UWISC  and  spatial  eigen\'alue 
results  fall  within  experimental  accuracy. 
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A  finite  element  method  has  recently  been  applied  to  the  steady  three-dimensional 
flow  in  a  spinning  and  nutating  cylinder*^  using  the  software  package  FIDAP.  FIDAP  is 
a  commerically  a^’ailable,  general  purpose  code  for  the  solution  of  incompressible  fluid 
flow  problems  governed  by  the  Navier- Stokes  equations.  As  is  the  case  for  finite-difference 
codes,  the  incompressible  continuity  equation  is  augmented  by  the  pressure,  which  has  a 
small  coefficient  10”®)  called  a  penalty  peu-ameter.  The  study  in  Reference  14,  used  ap¬ 
proximately  2200  elements  over  the  entire  cylinder  wdth  16  nodes  on  the  sidewall  boundary. 
In  the  general  finite  element  method  the  velocity  vector  is  approximated  on  each  element 
by  a  simple  polynomial  function,  but  in  this  case  the  velocities  were  approximated  by  lin¬ 
ear  interpolating  functions.  A  Galerkin  method  (or  weighted  residuals)  was  used  to  reduce 
the  Kavier-Stokes  and  continuity  equations,  together  with  the  boundary  conditions,  to  a 
system  of  nonlinear  algebraic  equations.  These  were  solved  by  a  quasi-Newton  iterative 
technique.  The  code  used  the  aeroballistic  reference  frame  (non-inertial)  for  steady  flow 
computations. 

Figure  3  shows  a  series  of  comparisons  between  UWISC  and  FIDAP  for  c/a  =  1.042,  a 
coning  angle  of  one  degree,  coning  frequencies  ranging  from  0.0385  to  0.0909,  and  Reynolds 
numbers  ranging  from  5  to  25.  FIDAP  data  were  obtained  from  Dr.  Simon  Rosenblat 
of  Fluid  Dynamics  International.  Proper  conversions  have  been  employed  to  express  all 
dimensionless  quantities  in  the  inertial  reference  frame.  For  Reynolds  numbers  ranging 
from  5  to  25.  both  the  UWISC  and  FIDAP  co'^es  show  consistency  between  roll  and  side 
moment  coefficients,  and  code  to  code  agreement  is  good.  A  constant  difference  of  about 
9.2%  in  the  results  of  UWISC  and  FID.\P  was  found  over  all  Reynolds  numbers  and  coning 
frequencies  considered. 

A  reasonable  range  of  dimensionless  parameters  must  be  prescribed  to  relate  the  com- 
ptited  outputs  for  code  to  codv  comparisons.  For  the  present  work,  the  physical  geometry 
is  restricted  to  that  of  a  completely  filled,  right  circular  cylinder.  Dimensional  analysis 
and  linear  theories  indicate  that  the  liquid  moment  coefficients  (roll  and  yaw  (normally 
called  side))  will  depend  upon  the  following  dimensionless  groups; 


Linear  Liquid  Momeni  Coejj.  -  F\Re,cJ a.T,kc/ 6]  (1) 

Given  the  use  of  a  SS-INS  code  that  retsuns  the  nonlinear  terms,  then  the  liquid 
moment  coefficient  will  also  depend  upon  the  precession  2mgle  =sinQc).  However,  a 
steady  state  code  would  require  —  0.  Hence,  the  present  codes  would  yield  a  dependence 
as  follows: 


SS  —  lyS  yorilinear  Liquid  Momeni  Coeff.  =  F[Re,c/a.T.  K^]  (2) 

The  case  for  low  Re  should  also  follow  this  formulation  and  will  be  examined  using  the 
dimensionless  groups  as  guides.  It  is  highly  possible  that  the  liquid  moment  coefficients 
are  linearly  related  to  q,-  for  <  20  degrees.  If  this  is  the  case,  then  two  of  the  remaining 

Rosentlai,  5..  Uooiir.g,  A.,  ar.d  EngUman,  M.  "rinite  £icmen*  Ca/cn/afton#  oj  VSicoe/ajitc  Flutd  Fiou  tn  c  Spinning 
(iTnf  AufaJinj  Ciixnde^.”  CRDEC-CR^d'^OiJ ,  Chemteal  RtiearcK,  Dctelopmeni  and  Engxneertng  Ctnur,  Aberdeen  Piovmg 
Greund.  M ar^iand.  December  J966. 
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three  parameters  can  be  held  constant,  while  the  behavior  of  the  liquid  moment  coefficient 
upon  the  third  parameter  can  be  explicitly  shown. 

Murphy  suggested  the  use  of  roll  and  side  moment  coefficients  for  small,  fixed 
precession  angles,  Uc,  defined  below: 

Roll  Moment  =  mLO^^'^lCiRMo  +  Tf^’c^LRM]  (3) 

Transverse  Moment  =  mia^^^T[CLSM  +  (4) 


where, 


mi 

a 

T 

Clrm 

ClRMc 

Cishf 

Cum 

Kr 

0c 


is  the  mass  of  liquid  in  a  fully-filled  container 
is  the  maximum  radius  of  the  container 
is  the  spin  rate  of  the  container  in  the  inertial  frame 
is  the  ratio  of  coning  rate  to  spin,  0c/0 

is  the  steady-state  liquid  roll  moment  coefficient  due  to  coning  motion 
is  the  liquid  roll  moment  coefficient  due  to  transient  liquid  spinup 
is  the  liquid  side  moment  coefficient 
is  the  liquid  in-plane  moment  coefficient 
is  sinoc,  where  is  the  precession  angle 
is  the  phase  angle  of  the  coning  motion 


Further.  Reference  11  gives  a  relationship  between  the  moment  coefficients  for  the 
linearired.  viscou'-  ’^•avier- Stoke;.-  equation^.  Hence,  for  small  precession  angles  and  all 
Reynolds  numbers, 

Clrm  =  -Clsm  (5) 

These  definitions  for  the  liquid  moment  coefficients  will  be  used  to  scale  computed 
lesuhb  and  cam  provide  for  comparisons  whth  various  sources  of  expevimf’nta)  data. 


III.  FINITE-DIFFERENCE  METHODS 


The  Navier-Stokes  equations  are  written  in  a  non-inertial  coordinate  system.  Descrip¬ 
tions  of  the  two  SS-INS  codes  (S.4ND  and  UWISC)  are  provided  in  References  8  and  9 
and  will  not  be  repeated  here  in  detail.  The  fundamental  approaches  of  the  two  codes 
are  similar,  however  the  solution  algorithms  are  different  (see  Section  II).  The  corriputed 
pressure  and  velocity  fields  are  input  to  sub-programs  that  compute  the  liquid-induced 
moments.  The  SAND  code  computes  the  pressures  and  shear  stresses  (using  a  second- 
order  finite-difference  on  velocity)  at  the  lateral  and  end  walls  and  then  integrates  these 
over  the  entire  cylinder.  The  UWISC  code  forms  an  expression  for  the  angular  momentum 
of  the  liquid.  This  expression  is  differentiated  with  respect  to  time  in  order  to  determine 
the  resultant  liquid  torques. 

The  SS-INS  codes  required  the  use  of  a  coordinate  system  where  steady  solutions 
would  exist.  Both  UWISC  and  S.AND  use  the  same  processing  coordinate  system.  Equa- 


tions  1  thru  5  are  defined  with  respect  to  the  inertial  spin  rate  tp.  Both  the  UWISC  and 
SAND  codes  used  components  of  the  intertial  spin  rate  to  describe  program  inputs.  The 
inertial  spin,  d,  is  the  sum  of  two  component  angular  velocities:  <p  =  (pp  +  (pcCusac  The 
UWISC  and  SAND  codes  define  the  dimensionless  variables  using  4>p  rather  than  d;  hence, 
the  notation  for  the  dimensionless  groups  must  be  augmented  with  a  subscript  p  for  the 
precessing  fratme  (7?ep,  Tp,  and  Cism^).  For  the  results  presented  in  this  report, 

proper  conversion  has  been  made  so  that  all  dimensionless  groups  are  defined  with  respect 
to  the  inertial  frame. 

Operational  differences  between  the  two  codes  will  be  discussed.  The  UWISC  code  was 
executed  on  the  Ballistic  Research  Laboratory  (BRL),  Launch  and  Flight  Di\'ision  (LFD) 
VAX  8600,  while  due  to  longer  CPU  times,  the  SAND  code  was  run  on  the  BRL  CRAY- 
XMP.  TABLE  1  documents  nominal  running  times  for  ’’converged”  solutions.  Although 
different  computing  machines  were  used,  run  times  are  expressed  in  VAX  CPU  units  witii 
1  hour  of  CRAV-XMP  time  equi\’alenced  to  10.6  hours  of  VAX  8600  CPU  units  (the  code 
was  not  altered  to  take  advantage  of  CRAY  vector  processing).  The  UWISC  code  uses  a 
typical  convergence  criterion  based  upon  the  Laplacian  of  the  solution,  while  the  S.\ND 
code  simply  recommends  that  additional  iterations  be  computed  to  assure  that  the  output 
reflects  a  ’'steady  state”  solution. 


Tablf  1:  Typical  8600  CPU  Times  (Hrs)  for  Solution  Convergence. 


Code  Grid  Points  Total  Reynolds  Number 


(r.O.z) 

Grid  Points 

10 

20 

30 

40 

UWISC 

11.12,42 

5,544 

0.738 

0.642 

0.621 

0.654 

SASD' 

11.24,21 

5,544 

2.100 

2.090 

2.090 

2.090 

*Keflects  appro.vin-iately  dO.OOO  iterations  per  solution  for  all  cases 


One  important  difference  between  the  codes  is  that  the  SAND  code  utilizes  dimen¬ 
sional  inputs,  whereas,  the  UY'ISC  code  calls  for  Re^.  c/a,  tmd  Qc-  A.  series  of  com¬ 
puter  runs  were  made  to  establish  dimensional  consistency  for  the  SAND  code.  In  these 
test  cases,  the  dimensional  variables  were  chosen  such  that  the  dimensionless  variables  did 
not  change.  Hence,  the  moment  coefficients  should  not  change.  The  results  &ie  shown  in 
T.A.BLE  2. 

The  computed  res..ltE  indicate  an  inconsistency  within  Clsm-  Clhm  involves  only  the 
viscous  stresses,  while  Cism  is  derived  from  both  viscous  ajid  pressure  terms.  Most  likely, 
the  computed  pressure  values  or  the  integration  of  these  pressures  along  the  cylinder  walls 
are  in  error  in  the  S.A.ND  code. 

Both  the  UWISC  and  S.UND  codes  have  been  extensively  used  to  compute  Clrsi  and 
Clsm  values  for  .6  <  Rt  <  45.  .03  <  r  <  .5,  and  1  <  c/g  <  5.2.  Figure  4  shows  a  typical 
comparison  of  Clrm  and  Clsm  between  the  codes.  The  Cism  values  for  SAND  have  not 
been  included  in  liglu  ot  tne  results  of  TABLE  2.  The  UWISC  code  is  consistent  for  Clrm 
and  Clsm  values.  Agreement  for  values  between  the  codes  is  shown. 


Table  2;  Dimensional  Consistency  V’erificalion  for  SAND 


(Hz) 

<Pc  (Hz) 

Viscosity  (cs) 

-Clrm 

Cisht 

Caee  A  {Rc  =  10,  c/a  =  4.32,  r  =  0.091,  Qe  = 

2‘) 

50.0 

4  55 

1.08x10® 

0.0237 

0.0417 

78.4 

7.13 

1.70x10® 

0.0237 

0.0391 

100.0 

9.10 

2.17x10* 

0.0237 

0.0387 

200.0 

18.20 

4.34x10® 

0.0237 

0.0310 

Case  D  {Rc  =  21.5,  c/a  =  1.042.  r  =  0.0797,  a. 

=  2’’) 

50.0 

3.99 

5.91x10^ 

0.0232 

0.0306 

100.0 

7.97 

1.18x10® 

0.0232 

0.0308 

200.0 

15.94 

2.36x10* 

0.0232 

0  0312 

Case  C  {Rc  ^  20.0,  c/a  —  5.20,  r  =  0.087,  = 

=  2“) 

50.0 

4  35 

4.15x10^ 

0.0296 

0.0422 

120.5 

10.48 

1.00x10® 

0.0296 

0.0498 

200.0 

17.40 

1.66x10® 

0.0296 

0  0505 

1.  CODE  LIMITATIONS. 

The  UWTSC  code  was  tested  over  a  wide  rajige  of  non-dimensional  coning  frequencies 
and  Reynolds  numbers  to  identify  code  limitations.  Figure  5  shows  the  results  for  a 
Reynolds  number  of  10,  cylinder  aspect  ratio  of  3  and  a  range  of  coning  frequencies.  The 
results  from  UWISC  are  compared  to  results  from  the  spatial  eigen\'alue  method  described 
in  section  I.  The  agreement  between  computational  methods  for  0  <  r  1  is  good. 
However,  the  UWISC  code  is  not  able  to  achieve  solutions  for  r  =  0  or  r  >  1.  For  r  =  0 
the  current  version  of  the  code  is  inoperable  since  the  Navier-Stokes  equations  have  been 
divided  by  r.  Solutions  generated  for  small  \-alues  of  r  near  zero  can  be  extrapolated  to 
T  =  0,  but  the  code  can  be  modified  to  accept  r  =  0  input.  For  t  >  1  the  code  would  not 
converge  for  a  wide  range  of  relcixation  and  convergence  acceleration  parameters.  It  is  not 
clear  that  this  limitation  in  the  numerical  scheme  can  be  corrected. 

Figure  G  shows  UWISC  results  for  a  cylinder  ^lspect  ratio  of  1.4SG,  .04  <  r  <  .1. 
and  values  of  Reynolds  number  to  300.  From  a  computational  standpoint  this  docs  not 
represent  a  Rc  limitation  in  the  UWISC  code.  Solutions  for  higher  Rt  can  be  achieved 
using  solutions  for  lower  Rc  as  initial  conditions.  However,  the  S.4ND  code  is  limited  to 
a  maximum  Rc  of  about  120  ®  and  this  limitation  may  be  a  result  of  the  explicit  solution 
algorithm.  The  spatial  eigen\’alue  method,  described  in  section  I,  has  been  run  for  these  Rc 
on  a  \'AX  SGOO  mini-computer  with  run  times  between  1  and  4  epu  minutes.  In  addition, 
the  spatial  eigenvalue  method  has  been  efficiently  run  for  Rc  as  large  as  3000,  with  no 
apparent  Re  limitation. 

The  real  limitation  in  achieving  high  Rc  solutions  using  UWISC  is  the  computational 
time  required  to  reach  a  steady  state  solution.  Figure  7  shows  the  average  run  time  per  r 
for  the  cases  plotted  in  the  previous  figure.  For  Re  =  300  the  run  time  on  a  CRAY  XMP/4S 
(without  using  vector  processing)  was  7  hrs.  In  this  case  the  solution  at  Re  —  200  was  used 
as  an  initial  guess.  Steady  state  results  for  larger  Re.  while  achierTible  using  the  solution 
algorithm,  may  be  prohibitive  in  light  of  the  computer  run  time  required  to  reach  them. 
In  this  regard,  the  sj^atial  eigenvalue  method  that  runs  very  efficiently  on  the  VAX  8600 


and  CRAY  computers,  would  be  the  preferred  choice. 

Several  additional  factors  should  also  be  considered.  The  results  of  Figure  6  were 
achieved  using  a  sirgle  computational  grid  for  which  the  grid  clustering  near  the  cylin¬ 
der  walls  was  increased  in  accordance  w’ith  the  wall  viscous  gradients  (i.e.  the  Reynolds 
number).  No  attempt  was  made  to  find  the  number  or  distribution  of  grid  pomts  that 
minimized  run  time  while  preserving  the  accuracy  of  the  solution.  In  tact,  the  number 
of  grid  points  used  in  each  case  was  chosen  to  exceed  that  required  for  a  unique  solution. 
A  strict  convergence  tolerance  («  =  1x10“®)  was  used  so  that  the  computational  results 
could  e%'aluate  the  predicted  relation  —Cirm  —  Cism,  to  within  l^'o  (on  a  VAX  8600  mini¬ 
computer  solutions  were  achieved  with  far  less  run  time  but  with  a  larger  value  for  f).  In 
addition,  fixed  \'alues  )f  the  relaxation  and  convergence  acceleration  parameters  were  used 
for  eacli  Rc:  there  was  no  attempt  to  find  the  most  efficient  choice  of  these  parameters.  It 
is  highly  possible  that  computer  run  time  for  a  single  case  (i?e,T,c/a)  could  be  reduced 
after  a  host  of  preliminary'  runs  had  been  made  to  optimize  these  factors.  In  practice,  the 
choice  of  relaxation  and  acceleration  parameters  are  made  based  upon  experience,  since 
the  uniqueness  of  the  solution  is  independent  of  these  choices.  The  computer  run  time, 
however,  is  dependent  on  the  particular  values  of  these  parameters. 

2.  NONLINEAR  EFFECTS  AT  LARGE  PRECESSION  ANGLES. 

T.A.BLE  3  shows  the  ratios  of  computed  veJues  of  the  moment  coefficients  for  precession 
angles  of  2  ami  20  degiees.  rt  oia  computed  data,  both  codes  indicate  that  nonlinear  effects 
for  Re  <  20  are  relatively  weak.  Hence,  the  use  of  a  linear  method,  such  as  the  spatial 
eigenvalue  approach,  can  be  used  with  confidence.  An  additional  interpretation  can  be 
made  with  respect  to  Equation  (5):  -Cmsf  cr  Cissf  for  precession  angles  less  than  20 
degrees  if  Rc  <  20. 


Tal.le  .3  Importance  of  Nonlinear  Effects 


Case 

Code 

CLs.\f  (20*^)/C/.s,tf  (2') 

Clrm  (20^)/CjffA(  (2'’) 

lif  -  20-0.  c/a  =  5.2U.  r  =  OsT 

I'WISC 

0.941 

1.000 

S.-\ND 

0.965 

1.004 

Rc  =  10.0.  c/a  =  -1.32.  r  =  .001 

CWISC 

0.941 

1.003 

SAND 

0.966 

1.004 

Rc  —  21.5,  c/a  —  I.-ISO.  r  =  .072 

UWISC 

0.941 

1.000 

SAND 

0.950 

1.000 

Rt  =  21.5.  c/a  -  1.0-12.  T  -  .079 

UWISC 

V.944 

1.004 

SAND 

0.958 

1.004 

IV.  CONCLUSIONS 


Two  incompressible,  steady  state,  fully  viscous  Navier-Stokes  codes  (SAND  and  UWISC 
have  been  compared  to  each  other.  The  following  conclusions  have  been  reached  to  date: 


1.  SAND  is  not  dimensionally  consistent  for  yaw  moment  predictions. 

2.  UWISC  and  SAND  are  consistent  for  roll  moment  predictions. 

3.  UWISC  verifies  that  for  small  Re  and  precession  angles  of  2“,  Clrm  —  ~Clsm- 

4.  For  coning  frequency,  UWISC  is  limited  to  |  t  |<  1 


5.  Without  optimizing  input  to  the  solution  dgorithm,  UWISC  code  computer  run  times 
are  prohibitive  for  steady  state  solutions  above  Re  =  300. 


6,  Linear  methods  can  be  confidently  and  efficiently  used  since  non-linear  effects  have 
been  shown  to  be  small. 
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Figure  1:  Compari.son  of  Ci/j.u  and  Cls^f  from  UWISC  code  with  Herbert’s  theory  (Rc  =  10). 
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Figure  4:  Comparison  of  Clrm  and  Clsm  from  VWISC  code  with  SAND  code  (c/a  =  4,32). 
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List  of  Symbols 


a  Maximum  radius  of  the  container 

c  Half-height  of  the  container 

c/a  Aspect  ratio  of  the  container 
Cum  Liquid  in-plane  moment  coefficient  (see  Eq.  4) 

Clrm  Liquid  roll  moment  coefficient  (see  Eq.  3)  using  inertial  spin  rate  {4>) 
CtRMp  Liquid  roll  moment  coefficient  using  non-inertial  spin  rate  {4>p) 

Clrmo  Liquid  roll  moment  coefficient  due  solely  to  spin  (see  Eq.  3) 

Clsm  Liquid  side  moment  coefficient  (see  Eq.  3)  using  inertial  spin  rate  (^) 
CisMp  Liquid  side  moment  coefficient  using  non-inertial  spin  rate  ((^p) 

Cp  Nondimensional  pressure  coefficient,  Cp  =  — 

Kr  sin(Qc) 

rriL  Mass  of  the  liquid  fill 

r  Radial  direction  normalized  by  a 

Re  Reynolds  number  {<pa'^/u)  using  inertial  spin  rate 

Rtp  Reynolds  number  using  non-inertial  spin  rate 

2  axial  direction  normalized  by  a 

Oc  Precession  angle 

€  Convergence  tolerance  for  UWISC  code 

w  Kinematic  viscosity  of  the  liquid  fill 

p  Density  of  the  liquid  fill 

T  Ratio  of  coning  rate  to  inertial  spin  rate  {<Pci<>\ 

Tp  Ratio  of  coning  rate  to  non-inertial  spin  rate  (^c/<^p) 

4>  Circumferential  angle 

cb  Spin  rate  of  the  container  in  the  inertial  frame  {<pp  +  (peCosOc) 

rpc  Coning  rate  of  the  ctintainer 

<pp  Eulerian  spin  rate 
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